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SUMMARY 

The computer simulation of thin film nucleatlon and growth, which was 
previously performed for the case of single monolayer, was modified to Include 
multilayer growth via the Volmer-Weber mechanism. The main conclusions drawn 
from the simulation results are: 

1. The kinetics of multilayer film growth Is nearly Identical to that of 
monolayer growth. 

2. When no reevaporation takes place, the cluster density resulting from 
multilayer growth Is higher at high coverage than that resulting from monolayer 
growth. 


3. When reevaporation does take place, the cluster density resulting from 
multilayer growth Is nearly Identical to that resulting from monolayer growth. 
This Is not due, however, to similarity In microstructure. 


INTRODUCTION 

Most studies of thin film nucleatlon and growth by means of computer sim- 
ulation reported In the literature (refs. 1 to 3) are limited to a single mono- 
layer, l.e., to the early stages of growth by the Frank-van der Merwe and 
Stranskl-Krastanov mechanisms. In these two mechanisms a monolayer of the 
deposited material Is formed before multilayer growth takes place. Only one 
study (ref. 4) which deals with multilayer growth has been published, and no 
comparison of the growth behavior under these two different modes of growth has 
been reported. 

Multilayer growth Is, however, very common and has been observed In many 
cases (see for example ref. 5). In fact, It constitutes one of the three basic 
modes of film growth, namely the Volmer-Weber mechanism, whereby multilayer 
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growth takes place from the very beginning. Thus, a study of the growth behav- 
ior via this mechanism, and especially comparison with the other mechanisms, 
seems to be warranted. 
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DESCRIPTION OF THE MODEL 

In order to allow comparison between the results obtained earlier (ref. 3) 
for the case of monolayer growth and those of the present study, the model 
employed here was made as similar as possible to the one employed there. Thus, 
the sample lattice selected for this study was a square lattice, and most of 
the basic assumptions are Identical, namely: 

1. Only single atoms Impinge on the surface 

2. Only single atoms (l.e., atoms which do not have any nearest neighbors) 
are capable of migrating on the surface. 

3. Mobile (single) atoms migrate on the surface by means of Jumps to 
nearest-neighbor sites and the Jump frequency Is constant. 

4. A cluster Is formed by two or more atoms occupying adjacent lattice 
sites. 

5. Only nearest-neighbor Interactions are taken Into account. 

6. No decomposition of clusters takes place. 

7. Only single atoms can reevaporate from the surface. 

8. The substrate surface Is free from sites of preferred nucleatlon, e.g., 
defects or Impurities. 

9. The boundary conditions taken are of an Infinite potential barrier, 

1 e., atoms cannot cross the boundary. 

For a discussion of these assumptions see reference 3. 


PROCEDURE 

lhe simulations were carried out on a microcomputer. A sample lattice of 
400 (20 by 20) points was used. Some test runs were also done on a 1600 (40 
by 40) point lattice, In order to check the effect of sample size. The results 
obtained In these test runs fell within the error range of the results obtained 
with the 400 point lattice. 

The x and y coordinates of the sites at which atoms Impinge on the 
substrate surface were chosen by means of a random number generator with homo- 
geneous distribution. An atom can Jump In four possible directions on a two- 
dimensional square lattice. The actual direction of each Jump was also chosen 
by means of a random number generator. In cases Involving reevaporation, the 
following procedure was used to determine whether a given (single) atom reevap- 
orates at a given time: a random number In the Interval (0,1) was generated. 
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If this number was smaller than the preassigned reevaporation probability, that 
atom was “erased." 

Eight simulation runs were carried out for each case studied, and the 
average as well as the standard deviation were calculated. In order to check 
the sufficiency of the number of runs, 16 simulation runs were performed for 
some cases. The results were within the error range obtained with eight runs 
and no reduction In standard deviation resulted. 


RESULTS AND DISCUSSION 

In order to study the effect of each of the processes Involved on the 
overall growth behavior, the simulation was done for each of the following 
cases. 


Impingement 

First, the case when the only process taking place Is Impingement of atoms 
on the surface was studied. This corresponds to the physical situation that 
exists when a film Is deposited on a very cold substrate. The three-dimensional 
microstructure resulting from this process on a 40 by 40 point lattice Is rep- 
resented In figure 1 for four values of the coverage. The percentage of atoms 
bound In clusters out of the total number of Impinged atoms Is given, as a 
function of coverage. In figure 2, where the results obtained for the monolayer 
growth are also given for comparison. It can be seen that the kinetics of 
growth In these two modes Is nearly Identical, although at coverage lower than 
0.4 a slightly higher fraction of the Impinged atoms Is bound In clusters In 
the case of multilayer growth. 

This Is not true, however, for the cluster density. From the results 
presented In figure 3, It Is seen that at coverages higher than 0.2 the cluster 
density resulting from the multilayer growth Is substantially higher than that 
resulting from monolayer growth. Thus, although these two mechanisms do not 
differ much In their kinetics, they do differ In the resulting mlcrostriicture. 


Impingement Combined with Surface Migration 

Next, the case when Impingement Is combined with surface migration was 
studied. This represents the situation when a film Is deposited on a substrate 
which Is hot enough to allow surface migration, yet Is too cold to allow reevap- 
oration of deposited atoms. The mlcrr structure resuming from this combined 
process on a 40 by 40 point lattice Is presented In figure 4. The kinetic 
results presented In figure 5 display a behavior similar to that obtained for 
the case of Impingement, namely the kinetics of growth In the multilayer mode 
Is almost Identical to that of the monolayer mode. Again, as In the case where 
only Impingement takes place, this Is not true for the cluster density. The 
results presented In figure 6 show that at a high coverage multilayer growth 
results In higher cluster density than that resulting from the monolayer mode. 
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Combined Impingement, Surface Migration, and Reevaporation 


Finally, the process of reevaporation of atoms from the substrate surface 
Is also Introduced. This represents the situation when a film Is deposited on 
a very hot substrate. For the sake of clarity, the results of only one such 
case, namely when the probability of each single atom reevaporating at any 
given time Is 0.5, are shown here. The mlcrostru.ture resulting from this 
process on a 40 by 40 point lattice Is presented In figure 7. The growth 
kinetics Is given In figure 8, and, as In the former two cases, It Is found to 
be almost Identical for the two modes of growth. 

Unlike the former two cases, however, the cluster density resulting from 
these two growth mechanisms, which Is given In figure 9, Is also similar. This 
similarity In cluster density Is not due, however, to similarity In microstruc- 
ture, since, as can be seen In figure 7, a considerable fraction of the atoms 
are located In "multilayer" positions. 


CONCLUSIONS 

The main conclusions of this study are: 

1. The kinetics of film growth In the multilayer mode Is nearly Identical 
to that of the monolayer mode. 

2. When no reevaporation takes place, the cluster density resulting from 
multilayer growth Is substantially higher, at high coverages, than that result- 
ing from the monolayer mode of growth. 

3. When reevaporation does take place, the cluster density resulting from 
the two modes of growth are very similar. This Is not due, however, to simi- 
larity In microstructure. 
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Figure 1. - Microstructure resulting torn direct impingement. 






Figure 2. - Percent o ( atoms nucleited directly upon 
impingement is i function of cover*, e. 



Figure 3. - Cluster density due to impingement is • function of coverige. 
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figure 4. Microstructure resulting from impinoement combined Aith surface migration. 
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Figure S. - Percent of nudeeted itoms upon impingement combined 
with surface migration. Impingement rite, 2. 5xl0 -3 atoms ( (site x 
time unit); Jump frequency, 1 per time unit. 



Figure 6. - Cluster density dye to impingement combined with surfice migration. 
Impingement rite, 2. 5xl0~ 3 atoms/isite x time unit); jump frequency, 1 per time 
unit 
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Figure 7. lAcroMructure resulting from combined impingement, surlace migration, and reevaporatcn. 





Figure 1 - Percent of nucleated atoms upon combined Impinge- 
ment, surface migration, and reevaporation. Implngemant 
rata. 2. 5xl0~^ atoms/tslte x time unit): jump frequency. 

1 per time unit; recvaporatlon probability, a 5. 



Figure 9. - Cluster density due to combined impingement, surface 
migration, and reevaporation. Impingement rate. 2. 5x10'* atoms/ 
(site i time unit); jump frequency, 1 per time unit; reevaporation 
probability. 0.5. 


